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MASS INFLUX OBTAINED FROM LLLTV 
OBSERVATIONS OF FAINT METEORS 

INTRODUCTION 


The present model of the meteoroid mass distribution is based on an 
extrapolation from ground based photographic observations of the larger 
meteoroids with masses of the order of grams to satellite-borne penetration 
measurements of meteoroids with masses in the microgram range. The mass 
range representing the greatest damage potential to manned space vehicles is 
from 1 to i 00 milligrams. The fact that the meteoroid population in this region 
must be inferred from an extrapolation over 6 orders of magnitude between 
two points that are determined by completely different properties of meteoroids 
through interactions that are poorly understood physically and cannot be ade- 
quately tested experimentally, has caused some concern among those respon- 
sible for establishing the meteoroid environment. Also, there are legitimate 
scientific reasons for extending the ground based optical measurements to 
fainter meteors. Of primary interest is the determination of the slope of the 
mass distribution curve or the population index parameter. This will greatly 
improve the confidence in the extrapolation as well as reduce the range over 
which the extrapolation must be carried out, and will provide a badly needed 
consistency check between the ground-based and satellite measurements. 

This paper discusses the techniques and results of using LLLTV obser- 
vations to determine the meteoroid mass distribution in the region from grams 
to milligrams. 


DESCRIPTION OF THE SYSTEM 


Recent developments in low-light-level television (LLLTV) systems 
have allowed the observation of much fainter meteors than could be photo- 
graphed. This improvement results primarily from the much higher quantum 
efficiency of the photodetector which results in much smaller integration 
times. This is particularly important in meteor work since it is advantageous 
to have the integration time shorter than the eventrduration. Our system con- 
sists of a Commercial Electronics camera chain using a Westinghouse WL- 
32000 Intensifier-SEC vidicon tube with a 105 mm f/. 75 Rayxar lens. This 
affords a 13° by 16° field of view. The effective integration time of the system 
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is very close to the standard frame time of 1/30 sec, which is ideal for 
meteor work. 

The ultimate theoretical sensitivity for the system is m^ = 14. 26 . 

This was estimated by requiring a star to produce 1 photoelectron at the photo- 
cathode per integration time. The system is invariably limited by sky back- 
ground which, even under ideal conditions, is 2 orders of magnitude above 

dark current. Stars as faint as m =11 have been observed which is close 

v 

to the theoretical limit of 11. 6 for a sky background of 300 m = 10 stars/ 
deg 2 , and a S/N = 5. 

For moving objects such as meteors, this limiting magnitude would 
apply to meteors moving nearly parallel to the optical axis so that they would 
remain within a resolving element for one integration time. For most meteors, 
the S/N will be decreased because the time they contribute to a resolving ele- 
ment is limited by the writing speed of their image. It is estimated that if the 
system has an observing limit of 11 magnitude for stars, it will see all meteors 
brighter than m = 6. 4 and 50 percent of the meteors brighter than m = 

8. 15. V V 


ANALYSIS PROCEDURES 


Even though the dynamic range of a TV image is limited, photometry of 
point images can still be performed over as much as 6 orders of magnitude by 
making use of the fact that the image spreads after it reaches saturation [ 1] . 
Thus, the amount of light associated with the image is a monotonic, if not 
linear, function of input. The difficulty lies in the fact that obtaining light 
curves from the TV monitor is a time consuming and laborious task, especially 
for the faint meteors. Until the special video processing systems presently 
being developed for this purpose are available, it will not be possible to obtain 
light curves on a sufficient sample of meteors to establish a good distribution. 

An alternative procedure, which is less time consuming, was adopted 
for interim use. This consists of treating the video system as a threshold 
detector and simply counting those meteors that are above the detection thresh- 
old. By varying the threshold through reduction of the lens aperture setting, 
a cumulative distribution in peak meteor magnitude is obtained. 

Figure i shows the results of this mode of operation during two observ- 
ing periods at Climax, Colorado. The camera was oriented toward the zenith. 
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Figure 1. Observed meteor rates as a function of limiting 
magnitude. (Error bars represent 1 c r limits based on the 
statistical sample. ) 





Care was taken to program the aperture settings to assure a uniform distribu- 
tion of observing time at each aperture setting throughout the night. The 
higher rates in the period from July 31 through August 9 may be attributed to 
the presence of the Perseids and the 6 - Aquarids during this period. The 
details of these data will be published separately. 

Taking the area of the sky within the field of view to be 640 km 2 the 
observed rates for the non-shower period can be expressed 

log <p = -15.482 + 0.5008 m v (l) 

where £ is the observed rate of meteors ( number /m 2 /sec) when m^ is the 

limiting magnitude. This limiting magnitude is taken to be the magnitude of 
the faintest AO star detectable on the video screen. In a sense, this limiting 
magnitude represents the faintest meteor that could be detected, although 
practically all detectable meteors will be brighter than this because of their 
writing speed. 


RELATIONSHIP BETWEEN PEAK LUMINOSITY AND METEOR MASS 


A statistical analysis by Jacchia, Verniani, and Biggs [2] of the peak 
brightness of meteors in terms of their mass, velocity, and entry angle yielded 
the result 

I = io- 4,63S m 0,9 v 3,5 (cos#) 0,6 (2) 

P o 

where I is the peak intensity in units of zero magnitude stars, m Q is the 

initial meteor mass in grams, v is the entry velocity in km/sec, and G is the 
entry angle or the angle between the velocity vector and the zenith. 

A similar result can be derived from the classical single body meteor 
theory by assuming that the rate of mass loss is equal to the energy input 
divided by the heat of vaporization L. In the free molecular regime 


. A pv 3 

m = 2L = 


S m ^ pv 3 


2Lpm 


% 


(3) 


where S is the shape factor (S = 1.208 for a sphere), p^ is the density of the 
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meteoroid, and p is the atmospheric density. This neglects radiation losses, 
which are small compared to the heat input even at the boiling point of Fe 
(3160° K) ; and neglects any change in shape factor with time. Neglecting decel 
eration, which amounts to only a small percent velocity change in the time it 
takes a small meteoroid to completely burn up, and assuming an exponential 
atmosphere, equation (3) becomes 


m 



2Lp 



'H 


-x/h 

e 


(4) 


where x = -vt cos 9 , h is the scale height, and p is the atmospheric den- 

sity at x = 0. This differential equation may be solved by separation of vari- 
ables. It is convenient to define t = 0 at meteor burnout, or m = 0. The 
solution is then 


m = 



(5) 


where 


n _ v cos 9 
P ~ h 

and m , the initial mass, is 
o 


Sv 3 p 


m =| 

o 


H 




6L p m ''0 
Differentiating the solution, 
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m=-30m o (l . (6) 

The peak m is obtained by equating m from (6 ) to 0 to find the time t^ 
when m is maximized. This yields 

p = 1/3 
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Putting this in (5), 

4 o 

m = - — m 
p 9 o 

The radiant intensity from a meteor is given by 


T T • 2 

I = - — mv^ 
2 


where r is the luminous efficiency. Using equation (7), 


. m v 3 

T T 4 . o - 

L = — - — - — cos 9 
P 2 9 h 


( 7 ) 


( 8 ) 


(9) 


The luminous efficiency for Fe has been determined experimentally from 
Trailblazer [3] and is expressed by t = 10~ 17 * 95 v (cm/sec). Figure 2 
compares the light curve obtained with this model with an observed light curve 
from one of the Fe Trailblazer meteors. 

For stony meteors. Cook, Jacchia, and McCrosky [4] recommend 

r = I0~ 18,91 v (cm/sec) (lo) 

which is consistent with estimates of Ayers, McCrosky, and Shao based on the 
Trailblazer measurements. 

Putting this value in equation (9), and choosing h = 5.4 km, which 
corresponds to 80 km, results in 

T = 10- 5,297 m v 4 cos 9 (ll) 

P o 

with the units the same as in (2). For the case of a 1-gm meteor at 
22 km/sec, typical of cases from which equation (2) was obtained, 

_ j 1.154 (equation 2) 

P ( 1.182 (equation 11) 

The fact that the simple theory yields almost identical results as the empirical 
approach, together with the desirability of having some theoretical basis for 
determining the functional relationship of luminous intensity with mass, 
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Figure 2. Comparison of theoretical light curve compared with the 
measured light curve of an artificial meteor. 



velocity, and entry angle, are the bases for choosing equation (ll) as the 
functional observing relationship for the analysis. 


SYSTEM RESPONSE TO MOVING TARGETS 


As was stated previously, the detection threshold of a moving object is 
increased because the photons are spread over a number of resolving elements 
instead of contributing to the signal in a single element. A first order attempt 
to derive the system response to moving targets is to simply require that the 
minimum detectable moving object deliver the same number of photons during 
the time it resides in a resolving element as a minimum detectable stationary 
object in one integration time. The residence time is 


t = 



r 

F co 


( 12 ) 


where A^ is the area of a resolving element, F is the focal length of the lens, 
and co is the angular rate. 

The criterion for detectability is 


I 


L 



r 


F co 


j t t 


(13) 


where I is the threshold intensity for stationary objects, I is the limiting 
T L 

intensity for moving objects, and t is the integration time. The photocathode 

is 40 mm in diameter, or 24 mm from the top to the bottom raster line. 

Since there are 525 lines, a resolving element is taken to be a square 24/525 

mm per side. Taking F = 105 mm and r = l/30 sec, 

T = I 76. 57 co ; co 2= 0. 01306 rad/sec (14) 

1j t 


Figure 3 compares this model of the response to moving point sources with 
measurements made in our laboratory. The measurements consisted of pro- 
jecting point sources onto a mirror mounted on a rate table which reflected the 
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Figure 3. Response of the SEC Vidicon to moving point images. (The 
solid curves represent the predicted response based on the 
detection criteria used in the analysis^) 
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image onto a projection screen. The point source consisted of a precision 50 n 
pinhole mounted in a 35 mm slide projector. Neutral density filters were used 
to vary the image intensity, and the various images were calibrated by a photo- 
multiplier photometer. The background lighting was provided by a small 
flood lamp with a controllable aperture. The angular rate of the mirror was 
adjusted for each filter until the image could no longer be detected on the 
monitor screen. 


The agreement between the simple model and the observed results 
justify the model. The only discrepancy occurs in the transition at the point 
where the image spends t in a resolving element. Since the image in the 
experiment is smaller than the resolving element, but is not an infinitesimal, 
some rounding off of the theoretical model in this transition region is to be 
expected, which of course is the observed result. 


Since the angular rate 

v sin 0 
u> = 


(15) 


and using r - 80 km typical of meteor heights, equation (15) becomes 


T = 


I 0.957v sin 0 


( 16 ) 


whichever is greater. 


RELATION BETWEEN INCIDENT FLUX AND OBSERVED FLUX 


The number of meteors observed per unit area time is given by 

oo co 

<p = f dfi COS0 / dv / \ dmn mv ( 17 ) 

hemisphere 0 m^, ^1 , 0 , y j 

where n is the directional mass velocity distribution (no. per unit area, 
mv 

time, solid angle with masses between m and m+dm and velocities between v 
and v+dv). The integration is carried out from the threshold m^ required to 
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produce an observed signal, which is a function of the threshold response of 
the system 1^,0, and v. Given only the observed 0 , clearly there is not 

sufficient information to solve the integral equation. Several simplifying 
assumptions are in order. First, to a good approximation the velocity and 
mass distributions are independent. Second, it will be assumed that the veloc- 
ities are isotropically distributed; i. e. , n is independent of 9 and 0 . 

mv 

Equation (17) becomes 


i r/2 

0 = 2tt J cos 9 
0 


sin 9 d0 f n dv f 


m T (i T 


e 


v) 


n dm 
m 


(18) 


Even assuming n^ is known, there is still insufficient information to define 

n . However, the fact that the observed 0 can be expressed as a power 
m 

law, equation (l), suggests the cumulative mass flux (no. per unit area 

time with mass m or greater) can be expressed 


N 


m 


7T f N dm = C 
J m 


-a 






where a is the population index. With this assumption equation (18) becomes 
n/2 00 _ v 

0 = 2C / cos 0 sin 0 d0 J n dvm^ “Lfi.v . (19) 

0 0 

The threshold mass is from equation (ll) 

m = I 10 5 * 297 v- 4 (cose)" 1 
T Xj 

But from equation (16) 


rn = 
T 


u 


T 

[ T 


1q 5. 297 y - 4 

(cos e)- 1 

o 

CL) 

M 

CD 

. * 

10 5.297 y— 3 

tan 9 

o 

VI 

CD 


( 20 ) 
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where 


0 = sin -1 — - — — 

o 0.957v 


Equation (19) becomes 


f = 2CI t -“ 10 - 5 ’ 297 “ 


CO d o ^ 

f v 4<X n dv f (cos0) 1+a sin0 d0 
L 0 V 0 


® ir/2 , , 

+ (0.957) -Q! f v 3a n dv f (cos 0) +a (sin0) ~° i d0 

0 V 0 (v) 

o 


( 21 ) 


It is convenient to introduce an average m defined as the mass of a just detect- 
able meteor having average velocity and 0 = 45° . From equation (20) 


m = 0.957 I T 10 5,297 v -3 tan 45° 


Equation (2l) becomes 


<P = 


2Cm ~ a 
3a 


CO 9 (v) 

(0.957) 01 J v^ 01 n dv f cos 0^ 01 sin 0 d0 
0 v 0 


CO 


7T/2 


+ J v 3ff n dv j (COS0) 1 01 (sin0) 
0 V 0 q ( v) 


1-Qf 


d0 


( 22 ) 


(23) 


The first integral is the contribution from those meteors moving nearly 
along the line of sight that remain in a single resolving element for one integra- 
tion time. Since 0 q is typically 3° , this contribution is small and can be 

ignored. The integral over 0 in the second integral must be evaluated 
numerically unless a is an integer. The lower limit 0 Q (v) is a function 

of v, as may be seen in Figure 4. This dependence is not strong for a » 1. 

Therefore an average value of v = 20 km/sec will be used which yields 

0 =3°. 

o 
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Figure 4. The value of the integral 1(0) in equation (25) 
as a function of 0 o and a. 






The velocity distribution was adopted from the work of Dohnanyi [ 5] , 
and expressed as 


j C N V ‘ 
11 = < 
v Ul. 


; 11.2 s -v ^ 16.6 


61X10 7 v- 4 ‘ 3 ; 16.6 < v ^ 72.2 


The normalization constant is C = 0. 001153. The (v) using this distribution 

N /x\ x 

is 19. 2 km/sec. The weighted average \y // (v) is shown in Figure 5. 
Equation (23) can be written 
/ 3ck 

<P = 2 1(e) N— ( 25 ) 

<v> 3 “ m T 

where N— is the cumulative isotropic mass flux of meteroids having mass 
m_ 

_ T __ 

m rp or greater, m is given by equation (22) , and 1(0) is 


l(0) = f (cos0) 1+a (sind)^~ a d0 
0 


-0.4m 


Since 1,^ = 10 


and the observational results were found to be 


(equation (l)] 


log <p = -15.352 + 0.5053 m 


, , . _ 0.5053 . T 

log $ = -15.352 - — ■ logI T 


Differentiating the log of equation (21) 


d (log <P) 
d (log'll) 
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From equation (26) a is found to be 1.263. Using this value in Figures 4 
and 5 


= 2.48 

1(0 ) = 0.933 

Combining equations (22), (25), and (26) results in 

log N — = -15.352 - 1.263 (log - 1.466) 

m T 

T 

-log 2 - log 2.48 - log 0.933 
or 

log N— = -14.24 - 1.263 log m_ 
ni T 

T 

The observed range from m^ = 7 to = 10 corresponds to a mass range 

of io -1, 334 to lo -2 * 534 gm. The result of equation (27) is compared with the 
existing distribution of meteors in Figure 6. The result from the July 31 - 
August 9 expedition is also shown, which contains the Perseids and 5-Aquarids. 
No attempt was made to alter the distributions in velocity and angle to account 
for the shower component. 

COMPARISON WITH OTHER DATA 

Figure 6 shows the relationship of the data obtained in this work with 
the current meteoroid mass distribution adopted by NASA [6] . Also shown 
for comparison is the Hawkins and Upton datum point [7] based on photographic 
meteors, and the points obtained from the Pegasus [8] and Explorer XXIII [9] 
penetration experiments. These data have been analyzed in terms of encoun- 
ter frequency, i. e. , number incident per unit area time on a surface without 
regard to angle of incidence of velocity. The NASA design curve refers to the 
number capable of penetrating a surface per unit area time, which can just be 
penetrated by a meteoroid with the specified mass under conditions of normal 
impact at the average velocity (assumed to be 20 km/sec) and average density 
(assumed to be 0. 5 gm/cm 3 ) . Integrating over the velocity and angular 


(27) 
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log N m (noJm z t*:) 


EXPLORER XXIII 
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Figure 6. Comparison of the results of this study with other work and with the adopted 
NASA meteoroid environmental design criteria. (The observed data represent 
encounter frequency, whereas the design criteria is weighted to express penetration frequency. ) 



distributions weighted appropriately for penetration mechanics results in the 
penetration frequency being less than the encounter frequency by a factor of 
approximately 2 [ 10] . 

The penetration data from Pegasus and Explorer XXIII were analyzed 
using recent calibration data [ 11] . The discrepancy between these data points 
and the NASA model is partially due to the difference between encounter fre- 
quency and penetration frequency, but also includes the departure from linear 
size scaling in thin metallic targets. Again, the design curve was derived 
from the actual penetration data using the conventional penetration formula, 
and will serve adequately as a design criteria so long as the same formula is 
used to convert back to penetration results. However, this departure should 
be considered in developing a true mass distribution. 


CONCLUSIONS 


A technique has been developed, using the SEC vidicon LLLTV system 
as a threshold detector for faint meteors, to obtain mass flux distribution data 
in the mass range from 1 to 100 milligrams. The analysis technique is based 
on peak intensities using the most recent values of luminous efficiency obtained 
from the Trailblazer measurements. The data are quite consistent with pres- 
ent photographic data at 1 gram, and the satellite data at 1 microgram, and 
tend to confirm the adopted NASA meteoroid model. 

It is recognized that these observations represent only one time during 
the year and may be subject to seasonal variations. Such effects are the object 
of a current investigation. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, December 10, 1971 
124-09-23-0000 
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